Feedback Effect on Landau-Zener-Stiickelberg Transitions in Magnetic Systems 
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We examine the effect of the dynamics of the internal 
magnetic field on the staircase magnetization curves observed 
in large-spin molecular magnets. We show that the size of the 
magnetization steps depends sensitively on the intermolecular 
interactions, even if these are very small compared to the 
intra-molecular couplings. 



Magnetization dynamics of nanoscale magnets, i.e. 
systems like Mni2-acetate and Fes have been studied ex- 
perimentally and theoretically lately . At sufficiently 
low temperatures quantum effects are observed, due to 
the discreteness of the energy levels involved. When the 
magnetization of a crystal of such molecules is measured 
during a sweep of the external magnetic field, a staircase 
hysteresis loop is obtained. The steep parts of the stair- 
case correspond to the values of the external magnetic 
field where there is a crossing of adiabatic energy levels. 
Several aspects of this quantum effect were studied in 
[p[-p^. In a zero-temperature calculation, one finds that 
the magnetization can only change in steps, very similar 
to the steps observed in recent experiments on high-spin 
molecules Mni2-acetate and Feg. At every crossing, only 
two levels play a role and the transition probability can 
be calculated using the Landau-Zener-Stiickelberg (LZS) 
mechanism |T^p^. Two parameters determine the LZS 
transition: The energy splitting at the crossing and the 
sweep rate of the magnetic field. 

The size of the energy splitting which leads to a LZS 
transition probability is determined by the off-diagonal 
terms in the Hamiltonian describing the system. A 
straightforward perturbative calculation shows that this 
splitting is roughly scales like r^'^™' where T determines 
the magnitude of the off-diagonal terms and Am de- 
notes the difference in magnetisation of the two relevant 
levels. In the absence of a transverse applied field the 
energy-level splittings in the high-spin molecules men- 
tioned above are so small that the probability for a sin- 
gle LZS transition is effectively zero, unless the applied 
longitudinal field is rather large (see for example [Q). 

In the crystal the magnetic field felt by a particular 
molecule is the sum of the external field and the internal 
field due to the presence of other magnetic molecules. As 



the inter-molecular magnetic couplings in these materi- 
als are weak compared to the intra-molecular interaction 
between the spins, it seems reasonable to consider the 
former as a perturbation. The purpose of this paper is 
to demonstrate that this argument fails in the case of 
LZS transitions. The point is that the LZS transition 
probability depends on the rate of change of the effective 
magnetic field at the crossing, which can be changed sig- 
nificantly by the presence of the internal magnetic field. 
The magnetization steps are found to be strongly affected 
by the type of interactions among molecules. We call 
this mechanism Feedback Effect on Magnetization Steps 
(FEMS). 

We first illustrate the effect for the case of the Mni2- 
acetate molecules. As a model Hamiltonian for this S = 
10 system we take R] 



Ti, = —Di S"^ — D4 (S^ + Sy 
—c t sin 6 Sx — {ct cos 9 



St) 



(1) 



Compared to the model of the extra feature in Hamil- 
tonian (|^) is the presence of a mean-field term, the 
strength of which we parametrize by A. It is clear that 
in this mean-field approach any new effect appears as a 
result of global changes of the internal field generated 
by all the molecules and is not due to local fiuctuations 
which should be treated separately ||l9t] . 

Quantitative results for the zero-temperature non- 
equilibrium dynamics of model can only be obtained 
through a numerical integration of the Schrodinger equa- 
tion. Using standard techniques we compute the 
magnetization steps for several values of A. The re- 
sults for Di = 0.64, Di = 0.004, tih angle 9 = 1° and 
sweep rate c = 0.001 (see |7j, we use dimensionless units 
throughout this paper) are shown in Fig. |l]. It is clear 
that the dynamics of the internal field can change the 
size of the magnetization steps considerably. FEMS is 
observed for all A 7^ 0. Note that the values of |A| we used 
are not unrealistic (|A| « D4 ^ Di), but rather small if 
we relate A to the dipole-dipole interaction which would 
yield a A which is 10 — 100 times larger p^ . 

At very low temperatures experiments^] show steps 
at lower values of H than the ones at which we observe 
steps in our calculation. In fact, for the set of model pa- 
rameters given in (|^) a much slower sweep rate c, much 
too slow for numerical calculations, is required if we want 
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FIG. 1. Magnetization dynamics of the Mni2-acetate 
model for several values for the intermolecular coupling, 
A = -0.005 (top curve), -0.003, -0.001, 0, 0.003 (bottom 
curve) . 
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-0.01 
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0.23 
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Step 2 


0.90 


0.78 


0.59 


0.48 


0.40 



TABLE I. 
steps in Fig. 



Transition probabilities corresponding to the 



to study the effect of the internal field at all level cross- 
ings. 

Therefore it is expedient to turn to a toy model in- 
spired by the one used to describe Fcg ||^. We take a 
5 = 2 model with the following Hamiltonian B: 



+TS,-ict + X{S,))S, 



(2) 



where we take D ^ 1, E ^ 0.08 and T = 0.08. These 
parameters are chosen such that we get two steps with a 
probability of about one half. 

In Fig. ^ we show the magnetization during a sweep of 
the magnetic field, with a sweep rate c = 0.01, for several 
values of A. We see that the FEMS effect is large. The 
transition probabilities are given in Table |. We clearly 
see a large change in the transition probabilities due to 
the presence of the internal field. 

A deeper understanding of the origin of the FEMS 
effect can be obtained by considering the system of N 
5 = 1/2 molecules described by the Hamiltonian 
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(3) 



where c is the sweep rate, T is the transverse field 
and J determines the interaction strength between the 



FIG. 2. Magnetization dynamics of the S — 2 model (|2|), 
for several values for the intermolecular coupling, A = —0.03 
(top curve), -0.02, -0.01, 0, 0.01 (bottom curve). 

molecules (|J| <C F). For simplicity we consider cou- 
plings between z-components only and assume the cou- 
pling between the molecules to be the same. Since |J| 
is small, we assume that we can make a mean-field-like 
approximation. The occurrence of FEMS does not de- 
pend on these simplifications (see below). This yields a 
Hamiltonian of a single molecule in a background field: 



n 



-T a-r 



{ct + X {(7z)) (7zj 



(4) 



where A oc J is an effective interaction. The system is 
prepared in the ground state, corresponding to a large 
negative time t, and the magnetic field is swept with 
constant velocity, until a large positive time is reached. 
Then, in the LZS case with A = 0, the transition prob- 
ability is given by the well-known LZS formula p = 
1 — exp(— TT F^/c). For A 7^ we write the Schrodinger 
equation corresponding to (^ in component form: 

iu'= (-c<-A(2|up-l)) u-Fd, (5) 
id' = (ci + A(2|up- 1)) d-Fu, (6) 

where we also have the normalization condition |up -I- 
|dp — 1. From numerical simulations we (see below) 
find that the tunneling is suppressed (enhanced) by the 
presence of a feedback term with positive (negative) A. 
This can be understood in terms of a changed effective 
sweep rate at the point of the transition. Because the 
effective magnetic field at the position of the molecule is 
given by ci-l- A (2 1), the effective sweep rate would 

be c -I- Xd{az)/dt. 

If A is small but non-zero, the mean field term only 
contributes at the point of the crossing. So we look at 
a Taylor expansion around the point of the transition tc 
(to be determined later). 



M(t) =Uo+Uiit- tc) + 0{{t - tcf ), 



(7) 



and a similar expression for d(t). We insert this expan- 
sion in and @ and obtain 



iu'{i) = -ciu{t) - Td{i), 
id'(t) = ci d(t) -Tu{i), 



(8) 
(9) 



with t = t — tc, c = c + AX ^(uq u*) is the renormalized 
sweep rate, and were 3?(u) is denoting the real part of u. 
We define tc as the point at which ct + X (cr^) changes 
sign, so 



tc^-{l-2\uo\ 
c 



(10) 



This enables us to write uq = uq and ui = ui. To deter- 
mine these constants we use Zener's solution fl^Jl^ and 
the properties of Weber functions. We find 
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-ttS/4: 



r(i + ^) 



(11) 



where p is the new probability for crossing, i.e. p = 
1 - exp(-7rrV£) and 6 = r^/£, and = A (1 - p)^ /c. 
The shift of the field at which the transition occurs can be 
written as AH = A (1 — p)^. To determine c we calculate 



71- A 

2 r(i + ii) r(i-2^) 

We find that (Oh can be approximated by 



5f(Mo Ml) 




(12) 



(13) 



with an error of maximally 10% (see Fig. ||). Within this 
approximation, c is given by the implicit equation 



c^c + XT^^/2^£e-^^/^'. 



(14) 



A simple relation can be obtained by replacing c by c on 
the right hand side. Then 



Ar^ y/2n/c 



„-7rr74< 



and 



p w 1 — exp 



c 1 + Ar2 y2V^e-r^V4< 



(15) 



(16) 



The resulting probabilities are shown in Fig. ^ The re- 
sulting probabilities based on a numerical solution of ([l^ ) 
or (|l^) for c show similar behavior. Also shown are the 
results obtained from the exact numerical solution of the 
Schrodinger equation (^. As a test of the validity of 
the mean-field approximation we also show the result 
of four interacting S = 1/2 spins, where we assumed 
A = (N — 1) J. Clearly the exact results confirm the 
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FIG. 3. Transition probability as a function of \/\/c, with 
r^/c such that in the LZS case, p = 1/2. The solid line is 
based on a numerical integration of the crosses are taken 
from a simulation of four interacting 5 = 1/2 spins and the 
dashed line is based on (p^. 



validity of the mean-field approximation and the simple 
analytic expression (|l6|). 

For values of A below approximately —2.0 the descrip- 
tion in terms of a renormalized sweep rate breaks down, 
which can also be seen from the singularity in the argu- 
ment of the exponent in eq. ([T^). This is because the 
picture of a simple, single crossing breaks down and the 
effective magnetic field at the position of the spin will no 
longer be a strictly increasing function of time. We con- 
clude that the expression ( p^ ) captures the main features 
of FEMS at a single crossing. 

The relevant parameters, controlling the size of FEMS, 
are F/^/c and A/^/c. Only for S=l/2 the energy-level 
splitting is directly proportional to F^ . For the high-spin 
molecules this is not case (see above), in particular for 
the levels with large |m|. Although in these cases the ef- 
fective energy level-splitting that enters the approximate 
two-level description can be small, a rather small value 
of A can nevertheless change the transition probability 
significantly. 

We have shown that the magnetization steps in the 
hysteresis loops of clusters of high-spin molecules may 
depend sensitively on the change of the internal mag- 
netic field at these steps. This implies that the dynamics 
of this internal fields has to be incorporated in a descrip- 
tion of the magnetization dynamics, even if its magnitude 
appears to be small compared to the other model param- 
eters (for large spin). At finite temperatures the effect 
described in this paper will be enhanced further due to 
the thermalization to states with lower energy and larger 
magnetization |20|. 
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